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Abstract 

We evaluate the invariant cross section for production of dileptons in forward rapidi- 
ties at RHIC and LHC, using the Color Glass Condensate formalism and present 
results for the nuclear modification factor Rd(p)A ^ ^ function of dilepton invariant 
mass for the most central deuteron (proton)-nucleus collisions. 



1 Introduction 



The recent results on suppression of hadron spectra in deuteron gold (dA) collisions in the 
forward rapidity region P of the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven 
National Lab (BNL) has generated a lot of interest in the applications of semi-classical QCD 
and the Color Glass Condensate 2J to RHIC. Even though the Color Glass Condensate 
is the prediction of QCD for the wave function of a hadron or nucleus at high energies 
m |3] , it is not a priori clear at what energy this happens. There is some experimental 
evidence jB] that RHIC may be at just high enough energy to see glimpses of the Color 
Glass Condensate. 

The applications of the Color Glass Condensate formalism to the heavy ion collisions 
at RHIC have been most successful at low pt IH] which probe the kinematic region where 
Xbj ~ 0.01 at mid rapidity. However, the mid rapidity region in heavy ion collisions is not 
the best place to look for the Color Glass Condensate because of the dominance of the 
final state effects, such as the energy loss of energetic partons JTj from the possibly formed 
Quark Gluon Plasma. 

The forward rapidity region in dA collisions (the deuteron fragmentation region) is the 
best place in a hadronic/nuclear collision to probe the Color Glass Condensate |B:. First, 
there is presumably no Quark Gluon Plasma formed in a deuteron gold collision so that the 
dominant final state effects such the jet energy loss from the plasma are absent. Second, 
the forward rapidity region probes the small Xbj part of the nuclear wave function and the 
large Xbj part of the deuteron wave function. This is the ideal situation for the Color Glass 
Condensate probes since the high gluon density effects in the nucleus which give rise to 
the Color Glass Condensate are the strongest in this kinematics. 

RHIC is a unique experiment in the sense that it has almost a continuous rapidity 
coverage, < y < 4, among its various detectors where it can detect various particles. 
The STAR detector can measure hadrons and photons in mid rapidity as well as at ?/ = 4. 
The PHENIX collaboration can measure hadrons at mid rapidity as well as dimuons in the 
rapidity region between 1.2 — 2.2 while BRAHMS has measured hadrons at mid rapidity as 
well as rapidities of 1, 2.2 and 3.2. Therefore, one has the chance to map out the rapidity 
{xbj) dependence of particle production and confront it with the predictions of the Color 
Glass Condensate formalism. Already, the qualitative agreement between the predictions 
of the Color Glass Condensate formalism [HI GUI ^] and the data from BRAHMS, both 
the suppression of RdA and its centrality dependence, are quite remarkable specially since 
all the available models in the market missed this suppression despite their many free 
parameters fTIl . 

Electromagnetic probes such as photons and dileptons provide another tool, in addition 
to hadrons, with which to investigate the properties of the Color Glass Condensate [T^ . 
They are cleaner in the sense that they do not undergo strong interactions with the other 
partons produced after the collision. Furthermore, in the case of hadrons, one needs to 
convolute the produced parton with the desired hadron fragmentation function. This means 
that a hadron measured at a given transverse momentum pt comes from the fragmentation 
of a parton at a yet higher transverse momentum kt = Pt/z. Since the collision energy is 
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fixed, this takes one to higher momenta which correspond to higher Xbj where the Color 
Glass Condensate effects become less dominant. Dileptons and photons do not suffer from 
this and are therefore, a better signature of the high gluon density effects. 

In this short paper, we provide a numerical analysis of the dilepton production cross 
section in deuteron-gold collisions at RHIC and proton-lead collisions at LHC in the 
forward rapidity region. We consider y = 2.2 at RHIC where the PHENIX detector is 
located and y = 5 at LHC where there are plans to measure dileptons. We show our 
results for the absolute cross sections as well as the nuclear modification factor Rd{p)A 
and show that, similar to hadrons, dilepton production is also suppressed in the forward 
rapidity region. 

2 Dilepton Production in d(p)A Collisions 

Our starting point is the dilepton production in quark-nucleus scattering using the Color 
Glass Condensate formalism as derived in [THj (see also [H] for an equivalent approach). 
The diagrams corresponding to this process are shown in Fig. where the virtual photon 
is emitted before, after or during the multiple scatterings of the quark from the target 
nucleus. 




Figure 1: Dilepton production in quark- nucleus scattering. 



It can be shown that the diagram where the virtual photon is emitted during the 
scattering of the quark from the nucleus is suppressed by the Lorentz factor 7 and does 
not contribute. Assuming the incoming quark is on shell and ignoring all the quark masses, 
one gets 
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[A;2 + M2(l - ^)] [(kt-zlty + M^{l-z)]\ j 

with It = qt+kt and qt and kt are the transverse momenta of the outgoing quark and dilepton 
pair respectively while z is the fraction of the energy of the incoming quark carried away 
by the virtual photon and M is the dilepton pair invariant mass, a is the Fourier transform 
of the dipole cross section given in (0). 

To proceed further, we will integrate over the transverse momentum of the dilepton 
pair. This can be done analytically and simplifies the expression for the cross section 
considerably. The transverse momentum integrated cross section is also more useful exper- 
imentally since the production rate is quite small, at least at RHIC. Rewriting the dipole 
cross section in the coordinate space, we get 
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To relate this to deuteron (proton)-nucleus scattering, we need to convolute Q with 
the quark (and anti-quark) distributions g(x, M^) in a deuteron (proton). As shown in 
[TI] . this can be written in terms of the deuteron (proton) structure function F2 
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and 



J2x[qf{x,M^)+qf{x,M^) 
f 



is the deuteron (proton) structure function with xp = '^[^^ ~^~^]- Note that the sum over 
quark and anti-quark flavors is different for protons and deuterons. Here, we will use eq. (0) 
to calculate the dilepton production cross section in deuteron (proton) nucleus collisions. 
In case of a proton projectile, we use the GRV98 parameterization [1^] of the structure 
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function . We note that, in this kinematic region, the incoming quark or anti-quark 
distributions in a proton are well known and there is very little difference between different 
parameterizations. For a deuteron projectile, we use the HKM parameterization 16 of the 
F2 structure function which includes shadowing. Again, since the projectile quarks and 
anti-quarks are in the large x^j region, the nuclear effects in the projectile deuteron are not 
large except at very large x^j ~ 0.7 — 0.9 where shadowing can be a 20% — 30% effect. 

The forward rapidity region in a d{p)A collision probes small x^j gluons in the nucleus. 
Therefore, the target nucleus is treated as a Color Glass Condensate. To proceed further, 
we need to know the dipole cross section <Jdipoie{xg, bt, rt). The dipole cross section satisfies 
the JIMWLK equation [T^. This equation has recently been solved on a lattice in [TH] . 
Alternatively, one can solve the large Nc limit of the JIMWLK equation, known as the BK 
equation ^H]- This has been done numerically by several authors [20\ . 

Alternatively, one can model the dipole cross section based on the known properties of 
the solution to the non-linear evolution equation in various limits. This has been done in 
[21] and a comparison to the HERA data on structure functions has been performed. It is 
shown in that one can fit the HERA data in the kinematic region < 50 GeV^ and 
Xbj < 0.01. Furthermore, this ansatz has a simple form which includes the right anomalous 
dimension in the extended scaling region [221 12H| , unlike some previous parameterizations of 
the dipole cross section [24 . Having the right anomalous dimension is extremely important 
specially since the recent data from BRAHMS at RHIC may indicate that we are in the 
region corresponding to the linear evolution and the BFKL anomalous dimension. 

In [21] parameterization, the dipole cross section has the following simple form 

J d\ CTdipoieixg, bt, n) = U (xg, nQs) (5) 

where 

M{xg,rtQs) = 1 - e-'^'"'^^''3= nQs > 2 

and 

nQs < 2 (6) 

The constants a, b are determined by matching the solutions at rtQs = 2 while 7^ = 0.63 
and K = 9.9 are determined from LO BFKL. The form of the saturation scale is taken 
to be = (xo/x)^ GeV^ with xq, A, Aq determined from fitting the HERA data on proton 
structure function F2. We refer the reader to ^T] for details of the fit. In case of a nucleus, 
we make the assumption that the saturation scale of the nucleus is = A^^^Q'^p. We 
now have all the ingredients necessary to evaluate the dilepton production cross section as 
given by eq. Q- 

It is instructive to consider the kinematic regions where a parameterization of the 
form used in [21] is appropriate. This parameterization is valid in the saturation region, 
defined by scales M < Qs, and in the extended scaling region defined by M < Qes where 
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Qes = Ql/Qo- In the case of the proton, the initial scale Qo is of the order of Aqcd while 
in the case of nuclei, it is Qo = Qs{xo) where Inl/xo is the initial rapidity where high 
gluon density effects become important. This is usually taken to be xq ~ 0.05 — 0.01. It 
is important to notice that the actual value of the extended scaling scale does not enter 
anywhere in our results, but that it just defines how high in dilepton invariant mass M we 
can go in this formalism. 

To get a feeling for the saturation and extended scaling scales, we note that at y = 2.2, 
the extended scaling scale of a proton is ~ 3 GeV for M ~ 3 GeV. This is why, at RHIC, 
we show results for dilepton invariant masses up to M ~ 3 GeV. For rapidity ?/ ~ 5 at 
LHC, the extended scaling scale goes up to 13 GeV. 

We show our results for the absolute cross sections, as given in eq. (jHI), in Fig. 0. For 
comparison, we show the cross sections for both proton-proton and deuteron-gold scattering 
for different dilepton invariant masses at rapidity y = 2.2, appropriate for the PHENIX 
detector at RHIC. 




Figure 2: Dilepton production at RHIC: y = 2.2 and bf = 0. 



It should be kept in mind that M ~ 3 GeV is quite likely at the edge of the extended 
scaling region for this kinematics while our formalism is best suited for the saturation 
and extended scaling region and could very well begin to break down as we go to higher 
dilepton masses. Also, since the saturation and extended scaling scales are larger for a 
nucleus than a proton, our results are probably more reliable for deuteron-gold collisions 
than proton-proton collisions for this particular kinematics. This will improve as we go to 
higher energies and/or higher rapidities such as those covered by LHC. 

The nuclear modification factor RdA defined as 

da'^^-^^^~^/dydM^d\ 
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is shown in Fig. Q for the most central deuteron-gold coUisions. We also show the nuclear 
modification factor for the proton-gold collisions for sake of comparison. It should be noted 
that nuclear shadowing of the incoming deuteron wave function is taken into account in the 
HKM parameterization ^H]- As is clear from the figure, there is a large difference between 
dilepton production in proton-gold and deuteron-gold collisions as there must be since the 
flavor and as a result, electric charge, decomposition of proton and deuteron are different. 
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Figure 3: RdA and RpA at RHIC: y = 2.2 and bt = 0. 

In Fig. (Q, we show the dilepton production cross section in proton-proton and proton- 
lead collisions at LHC, at the rapidity of y = 5. We now go to dilepton masses as high 
as M ~ 13 GeV because the saturation scale is now much larger so that our formalism is 
valid for higher dilepton invariant masses. 

In Fig. (jSJ, we show our results for the nuclear modiflcation factor RpA at LHC for the 
most central collisions proton-lead collisions. The suppression of the dilepton spectrum is 
much stronger than that at RHIC as expected since at LHC smaller Xf,j is probed. 



3 Discussion 

Semi-classical QCD extends the domain of applicability of weak coupling (perturbative) 
QCD to high energies where the naive perturbative QCD approach fails due to high parton 
densities. The resulting state of a hadron or nucleus at high energy where parton densities 
are high is called a Color Glass Condensate. While the existence of such a high parton 
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Figure 4: Dilepton production at LHC: y = 5 and bt 
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density state follows from QCD, the energy at which this happens can not be derived, at 
the moment, from the theory itself and needs to be determined experimentally. 

While there is some evidence in favor of Color Glass Condensate from HERA on 
electron-proton Deep Inelastic Scattering |21j 24j, the BRAHMS collaboration at RHIC 
may have the best signature of the Color Glass Condensate in a nuclear environment so far 
in their measurement of the negatively charged hadrons in the forward rapidity region. To 
verify that this indeed the case, it is important to investigate the predictions of the Color 
Glass Condensate formalism for other processes such as dilepton production. 

In this paper, we provide predictions for dilepton production cross sections in proton- 
proton and deuteron (proton)-nucleus collisions in forward rapidity regions at RHIC and 
LHC. The universal ingredient in this cross section, as well as cross sections for hadronic 
observables, is the quark anti-quark dipole- nucleus cross section which is subject to the 
non-linear evolution equation and which has recently been solved. In this work, we use an 
ansatz for the dipole cross section which has been successfully used to fit the HERA data 
on proton structure function. This model includes the physics of the BFKL anomalous 
dimension as well as scaling properties of the Color Glass Condensate and therefore, is well 
suited for our purpose. 

There are several caveats in this work which need a more careful treatment than con- 
sidered here. First, the dipole ansatz as given in (jH)) has been shown to work for a proton 
target, but has not been used for nuclei. In this work, we assumed that the A dependence 
comes in through the saturation scale Q^A^^^. Ideally, one would like to have this only in 
the initial shape of the dipole. The subsequent A dependence will then be determined by 
the non-linear equations. There are some indications that as long CLS Xfyj IS not too small, 
the assumed A^^^ dependence may be fine for large nuclei j^. 

Another caveat is that the dipole model used here does not have the correct high pt 
behavior. In other words, it does not match onto the Double Log DGLAP limit which is 
contained in the Color Glass Condensate. This may not be very important in the forward 
rapidity region since the pt coverage in this region is limited due to kinematics of the 
experiment. 

Most importantly, it remains to be seen whether dileptons at low invariant mass (below 
J/ip) peak can be measured at RHIC. This is most important since this is where our pre- 
dictions are most reliable at RHIC. If and when RHIC future upgrades allow measurement 
of dileptons at higher rapidities, one can go to higher dilepton masses since the saturation 
scale is larger at more forward rapidities. Therefore, it is highly desirable to have the 
capability to measure dileptons (as well as photons) at RHIC in as forward rapidity as 
possible. Also, the STAR collaboration may be able to measure direct photons as well as 
photon -|- jets at rapidity y = 4 in the next deuteron-gold run at RHIC. The photon + 
jet process is unique in the sense it directly probes the dipole cross section. Color Glass 
Condensate predictions for this process will be presented elsewhere [21] ■ Finally, having 
another deuteron-gold run at RHIC in the near future will establish conclusively whether 
the Color Glass Condensate has been observed at RHIC, which seems to be the case, and 
will allow one to investigate its properties in detail. 
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